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Abstract. We investigate a space-inhomogeneous discrete-time quantum walk in one dimension. We siiow that the 
walk exhibits localization by a path counting method. 
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1 Introduction 

As a quantum counterpart of the classical random walk, the quantum walk (QW) has recently attracted 
much attention for various fields. There are two types of QWs. One is the discrete-time walk and the other is 
the continuous-time one. The discrete-time QW in one dimension (ID) was intensively studied by Ambainis 
et al. [I]. One of the most striking properties of the ID QW is the spreading property of the walker. The 
standard deviation of the position grows linearly in time, quadratically faster than classical random walk. 
The review and book on QWs are Kempe Kendon [3], Venegas-Andraca [1], Konno [5], for examples. 

In the present paper we focus on discrete-time case. The model considered here is a space-inhomogeneous 
two-state ID QW. The two-state corresponds to left and right chiralitics defined in the next section. Let 
Pn{0) denote the probability that the walker returns to the origin at time n. The model is said to exhibit 
localization if lim„^oo P2n{0) > 0. The homogeneous two-state ID QW except a trivial case does not exhibit 
localization, see [1], for example. The decay order of p2n(0) is closely related to the recurrence. As for the 
recurrence property of QWs, see Stefanak et al. [H1[I1IH]- Localization of the homogeneous model was shown 
for a three-sate ID QW in [5], a four-state ID QW in [TU], and a multi-state QW on tree in [TT]. Mackay et 
al. [H] and Tregenna et al. [13] found numerically that a homogeneous 2D QW exhibits localization. Inui et 
al. [2] and Watabe et al. [TS] showed the phenomenon. In higher dimensions, a c?-dimensional homogeneous 
tensor-product coin model does not exhibit localization [7] . Oka et al. 1 16 analyzed localization of a two-state 
QW on a semi-infinite ID lattice, which is closely related to the Landau-Zener transition dynamics. Through 
numerical simulations, Buerschaper and Burnett [T7] and Wojcik et al. [TH] reported that the dynamics of 
the two-state ID QWs exhibits from dynamical localization, spreading more slowly than in the classical 
case, to linear diffusion like the homogeneous two-state ID QW as the period of the perturbation is varied. 
Linden and Sharam |19) investigated a similar inhomogeneous two-state ID QW where the inhomogeneity 
is periodic in position. They showed that, depending on the period 2k, the QW can be bounded for even 
k and unbounded for odd k in time. The former case corresponds to localization and the latter case to 
delocalization. An interesting question is whether localization emerges even for a simpler inhomogeneous 
two-state ID QW compared with the previous models. The present paper gives an affirmative answer to the 
question. Our result could be useful for quantum information processing by controlling the spreading of the 
walker. 

The rest of the paper is organized as follows. Section 2 gives the definition of our model. In Sect. 3, we 
present our main result (Theorem 13.21) of this paper. Section 4 is devoted to the proof of Proposition 13.11 
In Sect. 5, we prove Theorem 13.21 Finally, we consider some related models in Sect. 6. In contrast to our 
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model, we show that a similar simple inhomogeneous two-state ID QW, whose limit theorem presented in 
[20] . does not exhibit localization. Moreover, for the corresponding classical random walk also, localization 
does not emerge. 



2 Definition of the walk 



In this section, we give the definition of the inhomogeneous two-state QW on Z considered here, where Z is 
the set of integers. The discrete-time QW is a quantum version of the classical random walk with additional 
degree of freedom called chirality. The chirality takes values left and right, and it means the direction of the 
motion of the walker. At each time step, if the walker has the left chirality, it moves one step to the left, 
and if it has the right chirality, it moves one step to the right. Let define 



\R) 



where L and R refer to the left and right chirality state, respectively. 

For the general setting, the time evolution of the walk is determined by a sequence of 2 x 2 unitary 
matrices, {Ux ■ x e Z}, where 



Ux 



ax bx 
Cx dx 



with ax,bx,Cx,dx G C and C is the set of complex numbers. The subscript x indicates the location. The 
matrices Ux rotate the chirality before the displacement, which defines the dynamics of the walk. To describe 
the evolution of our model, we divide Ux into two matrices: 



Px 



with Ux = Px + Qx- The important point is that Px (resp. Qx) represents that the walker moves to the left 
(resp. right) at position x at each time step. 

For a given sequence {iOx : x € Z} with Wx G [0, 27r), our previous paper [20] treated the following Ux'- 







■ " 





■) Qx — 





Ux = Ux[i^x) = 



1 



1 -er'^ 

In the present paper, for a given sequence {lJx ; x G Z}, we consider 

1 



Ux = Ux{uJx) 



V2 



1 

-iujj- 



-1 



(2.1) 



(2.2) 



In particular, here we concentrate on a simple inhomogeneous model depending only on a one-parameter 
uj e [0, 2tt) as follows: 



f/o = ?7o(c^), Ux = Ux{Q) iix^Q. 



(2.3) 



So when a; 7^ 0, our model is homogeneous except the origin. If = 0, then this model becomes homogeneous 
and is equivalent to the Hadamard walk determined by the Hadamard gate Ux = Ux{0) = H: 



H 



1 
71 



1 1 
1 -1 



In this paper, we take (p^, — "^[1/V2, */\/2] as the initial qubit state, where T is the transposed operator. 
Then the probability distribution of the Hadamard walk starting from tp^, at the origin is symmetric. 

Let S„(Z,m) denote the sum of all paths starting from the origin in the trajectory consisting of / steps 
left and m steps right at time n with n — I + m. For example, 

S2(i,i) = gPo + PQo, 

S4(2, 2) = Q^PPo + P^QQo + QPoQPo + PQoPQo + PQoQPo + QPqPQo- 
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The probability that our quantum walker is in position x at time n starting from the origin with 
^[l/%/2,i/%/2]) is defined by 



|2 



P{Xn = x) = ||^„(Z,m)(p* 
where n = Z + to and x = —I -\- to. The following is the important quantity of this paper. 

p„(0) = P(X„ =0). 

This is the return probability at time n. Remark that p2n+i(0) = for n > 0. For our model with w = tt, a 
direct computation implies 

2 10 40 

P2(0)-^-0.5, p4(0) = ^ = 0.625, p6(0) = - = 0.625, 

P8(0) = ^ = 0.66406..., pio(O) = ^ = 0.66406..., pi2(0) = ^ = 0.63476 ... . 

In fact, as a consequence of our main result (Theorem 13. 2p . we have lim„^.oo P2n(0) = (4/5)^ = 0.64. 
Therefore the QW with uj = n exhibits localization. On the other hand, for the Hadamard walk case (i.e., 
iu = 0), 

pf^(O) = 1=0.5, pf)(0) = 1=0.125, pf)(0) = A^O.125, 

Pi"^ (0) = P = 0.07031 . . . , (0) = 1^ = 0.07031 . . . , p^f ^ («) = = 0-04882 .... 

Superscript {H) denotes the Hadamard walk. In this case, it is known that lim„_^cio P2n(0) = 0, (see Sect. 
6). So the QW with w = does not exhibit localization. 



3 Our result 

In this section, we present our main result on the inhomogeneous two-state ID QW. Let 



*2„(0) 



(3.4) 



for n > 0. This is the probability amplitude at the origin at time 2n, where the upper (lower) component 
corresponds to the left (right) chirality. Remark that «'2„+i(0) = ^[^'2^+1(0), ^-^^^1(0)] = ^[0,0] for n>Q. 
Let P = {1, 2, . . .}. Then we have 



Proposition 3.1 For n>l, 



k=l 



(ai,...,afc)eP^ 
a\-\ (-a/j— n 



where 



7± = 
C± = 



2,. ,. , „:„2 



— coso; ± ivl + sm" w, /i±=smcj=p 
W 2 |(1 + sin^ io) =F sinw + sin^ Lu I , 



sm w. 



E 

n=l 



-i-z^ + vTT^ 



The proof is given in Sect. 4. By using this proposition, the following main result of this paper can be 
obtained. As for the proof of the theorem, see Sect. 5. 
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Theorem 3.2 For our inhomogeneous two-state ID QW with the parameter uj G [0, 27r) defined by Eqs. \2. 2\ 
and \2.3\. we have 

lim P2„ = ^— =: c , 

n^oo \ 3 — 2 COSCJ J 

where P2n(0) is the return probability at the origin at time 2n. 

We present some properties on the above limit c{ui) (see Fig. 1). (i) c(w) = c{2tt — uj). (ii) c{uj) is strictly 
increasing in w G [0,7r]. (iii) c(0) = < c{uj) < c(7r) = (4/5)^ for any w G [0,7r]. Therefore if the model 
is inhomogeneous, i.e., uj G (0,27r), then it exhibits localization, i.e., c{lu) > 0. When w is the uniform 
distribution on [0, 27r), we see that i?[c(-)] = (25 — 7-\/5)/25 = 0.3739 . . ., where £'[c(-)] is the expectation 
of c{uj). As we will discuss in the last section, for another inhomogeneous two-state ID QW defined by Eq. 
12.11 we have lim„„^oo P2n{0) — 0. That is, the QW does not exhibit localization. 

QG 
05 
04 
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Figure 1: The plot of c{uj) 



4 Proof of Proposition 13.11 

In this section, we prove Proposition 13 . II bv using a path counting approach. To do so, we first consider the 
Hadamard walk starting from location m (> 1) on Z_|_ = {0, 1,2,.. .} with an absorbing boundary at (see 
Ambainis et al. [T] for more details, for example). In this model, Pq and Qo do not appear, since we consider 
only {Ux = H : X > 1} . Therefore the definition of the walk yields 



V2 



for a; > 1. Then 



1 
3 



1 1 







1 -1 



(x> 1). 



Let be the sum over possible paths for which the particle first hits at time n starting from m. For 

example. 



^(CX),1) 

"5 



P^QPQ + P 



We introduce R and S as follows: 



R 



1 

71 



1 -1 





s 



1 

71 





1 1 



We should remark that P,Q,R and S form an orthonormal basis of the vector space of complex 2x2 matrices 
with respect to the trace inner product {A\B) — tr{A*B), where * means the adjoint operator. Therefore 
^(oo,m) -|_^^ written as 



4 



Noting the definition of we see that for m > 1, 

c(oo,m) _ '^{'>o,m-l) p ,^(cx3,m+l)^ 

Then we have 



„(oo,m) _ 1 (oo,m-l) 1 (oo,m-l) „(oo,m) „ _ 1 (oo,m+l) 1 (oo,Tn+l) 

~ /K t'n-l "I" /Tr'n-l ' 'in ~ /?; ^n-l /K n-1 i 



^(oo,m) _ (oo,m+l) \_ (oo,m+l) „(oo,m) _ J_ (oo,m-l) ^ (oo,m-l) 



From the definition of S^(°°''"'', it is easily shown that there exist only two types of paths, that is, P . . .P 
and P ...Q. Therefore we see that q^'"^^ = si°°'™^ = for n > 1. We introduce generating functions of 
Pn and r„ 'as follows: 

oo oo 
n=l 71=1 

Then we get 

p(oo,m)(^) ^ ^ ^(oc,m-l)(^) ^ ^ r (z) , 

V 2 V 2 

v 2 v 2 

Solving these, we see that both p(°°'™)(z) and r(°°'"H-2) satisfy the same recurrence: 

r(°°'™+2)(^) + \/2 - r(°°''"+i)(^) - r(°°'"'(2) = 0. 
From the characteristic equations with respect to the above recurrences, we have the same roots: 



A± = — 



V2z 

The definition of sl°°''^'' gives pl^'^"^ = (n > 2) and p^""'^-* — 1. So we have p^°°'"^^(2;) = z. Moreover noting 
liuim^oo P^°°'"^\z) < oo, the following explicit form can be obtained: 



Therefore for m = 1, 



-1 + vT" 
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Next we consider the Hadamard walk starting from location m{< —1) on Z_ = {0,-1,-2,...} with an 
absorbing boundary at 0. Let 2^"°°'™^ be the sum over possible paths for which the particle first hits at 
time n starting from m (< —1). Similarly we see that 

z 

So for m = — 1, 



1-VTV. 



4 
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Remark that ri°°'^^ + Sn = for n > 1. Let S+ = !E;i°°'^^(5o and S„ = °°' ^^Pq, where 



1 e"^ 








That is, ^+ (rcsp. ^„ ) is the sum of all paths for which the particle first hits at time n starting from the 
origin restricted in the region Z_|- (resp. Z_). Therefore we obtain 



Lemma 4.1 (i) If n>4, and n is even, then 



Et = ri^t^ RQo = 



„(oo,l) 



where 







n=l 



A-oo-1) 



SPo 



J-oo-l) 
'n-1 



-oo,-l)^n ^ 1 - VI + 





1 e''^ 



fwj Ifn is odd, then 







7T+ _ 





1 e*" 








Put S* = S+ + S„ . Prom this lemma and Sn °°' = — ri°°'^^for n > 1, we have 



1 



where 



(-1)^ 



if n = 4m — 1 and m > 1 



_i (2to-1)! 

22'"-i(m- l)!m! 

if n 7^ 4m — 1, n > 2 and m > 1 

-1 ifn = l 



In fact, 



r* = -1, r* = 0, r* = 1/2, = = r * = 0, = -1/8, r* = = = 0, . . . . 
Then the generating function of r* is as follows: 

n=l ^ 

The definition of S* yields 



fe=l (ai,...,afc)eP'=: \i=l 
aiH |-a(s=n 



*-(o)=E E in--. 

fe=l (ai 
ai 

where P = {1, 2, . . .}. Then a little algebra gives 

1 

■ \ — cl — T+ + — T-/ 



e"'" 1 




" 1 " 


-1 -e''^ 


^/2 





1 

V2 
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where 



7± = — cos 



sin^ u, /i± = sincj =p V 1 + sin^ cj, 



C± = a/2 |(1 + sin^ w) =F sinwVl + sin^w|. 



Noting that 



n^2a, I V'' = 



1 




1 




■ 1 " 


2fe 


-1 




V2 





we have the desired conchision. 



5 Proof of Theorem [33] 

In this section, we prove our main resuU, i.e., Theorem l3.2l By Proposition 13. 11 we will compute generating 
function of *i^^(0). Put 

— ^2n-i ^± — 7±/2. Then we see that 

oo 



'"""E^E E in-". 



1 



oo ( n 



a 



n=l U'=l (ai,...,afc)6P'=: V = l 
diH htifc— n 



C2 



OO ( n 



EE E n-". 

n=l U=l (ai,...,afc)eP''': \i=l 
aiH \-ak=n 



1 

71 



1 - M+ 



oo f oo 



EE E n-". 

fe=l L„=fc (ai,...,afc)6P'=: \J = 1 
aiH |-afc=n 



OO f OO 



c^EE E n-". 

fe=l Ui=fe (ai,...,afc)eP'=: \i=l 
aiH hafc=n 



1 

71 



1 - 

CI 



OO ^ ^ oo 

E {(-1 - + ' + E {(-1 - + V^)--} 



1 {-I - + ^/l + ^)u+ (-1 - z2 + 7T^^)u_ 



s/2\ Cl l-(-l-z2 + VTT^)?i+ C*! 1 - (-1 - z2 + ^/T~^?)^i- j ' 
The first equality comes from Proposition 13. II As for the third equality, we should remark that for fc = 2, 



E E 

"=2 (ai,a2)eP^: 

£11+02=1 



2n 



-1 - z2 + Vl + z4 



= (-1 - z2 + v/l + z4)2 
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In the similar fashion, for general A; > 1, we have 



■r^=^ (ai,...,afc)GP'': \i=l 
Q-i -| \-ak—n 

Noting that the initial state *^^^(0) = l/\/2 and 



C2 



C2 



we obtain 



n=0 



where Z = — 1 — + v^TT^- Next we consider generating function of ^'^^■'(O). From the initial state 
*(^)(0) = i/V2and 



CI 



Cl 



we similarly get 



Therefore we have 



2 + Z cos a; 



h h x/2(2 + 2Zcosu; + Z2)' 

" ^ ^ x/2(2 + 2Zcosw + Z2)' ^ " x/2(2 + 2Zcosa; + ^2) ' 



(1 — smco)Z 



where ^lf^'^\o) (resp. ^'i'^'^^O)) is the real (resp. imaginary) part of ^^\o) for A = L,R. A direct 
computation gives 

oo oo 

^*^^)(o)." = y*(«'^)(o)." 



n=0 



n=0 



4 — Scosw + 2(1 — cosa;)2z2 _(_ ^2 — cosa;)z^ + (2 - cosw)(l + z2)vTT^ 



yv!>^^)(0)." = - 



n=0 



2\/2 {3-2 coscj + 2(1 — coscj)2z2 + (3 — 2 cosa;)z^} 
(1 + sinw) {1 + 2(1 - cosw)z2 + + (-1 + z^)Vl + z^} 
2\/2 {3-2cosw + 2(l-cosw)202 + (3_2cosa;)24} ' 



^ " ^ 2\/2 {3- 2cosa; + 2(1- cosa;)2^;2 + (3 _ 2 008^)^4} 



n=0 

Then we obtain 



*rH0)=v,.r'(0)-4^^^^cos(n.0), 



3 — 2 cos w 



^/2(l - cos g;) (1 + sing;) 



*^„'"^(0) 



(3 — 2 cos Lu) vl + si 



2 

sm uj 



sin(n^o)j 



ir,(.RM)rn\ ^2(1 - cosw)(l - smw) . , 

*2n (0) ~ : / . o sm(n6lo), 

(3 — 2cosw)v 1 + sin oj 
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where sin 6*0 = (2 — cosw) vT+sii?aj/(3 cosw — 2), cos^o = —(1 — coscj)^/(3cosa; — 2) and /(n) ~ g{n) 
means f{n)/g{n) — > 1 as n ^ oo. Concerning the above derivation, see pp. 264-265 of [H], for example. The 
definition of P2n(0) gives 

so the proof of Theorem 13.21 is complete. 



6 Discussion 

In the last section, we consider some relations between our model and other related ones. For the Hadamard 
walk (homogeneous model), a similar argument yields 



oo oo 

h h 2V2 V vT+z^ 



-l + z 



Ti— n— V \ VI 



Therefore we have 



irn 

As for the result, see [T], for example. Then lim„_j.ooP2rf''(0) — 0- Moreover, using Proposition 4.3 of [20], 
we obtain P2n(0) = P2n v^)- Here P2n(0) is the return probability at the origin at time 2n for another 
inhomogeneous model defined by Eq. 12.11 which was studied in [20]. Therefore the QW has also the same 
decay order as the homogeneous walk, i.e., Hadamard walk: 

pUo) - — • 

nn 

So lim„_^oo P2n(0) = 0, then the QW does not exhibit localization. This is in great contrast to our model. 
For the inhomogeneous classical random walk starting from the origin on Z, we similarly get 



p q 



(6.5) 



where pl^"* (0) is the return probability at time n for the classical walk. In this model, a walker at location x 
moves one step to the left with probability and one step to the right with probability Qx where Px + qx ^ ^ 
for any x S Z and Px = p, qx = q for x G Z \ {0}. From Eq. 16.51 we have 



Therefore lim„_^ooP2ri (0) ~ ^^^d localization does not occur. If p ^ q, thenp2n(0) decays exponentially. If 
p = q=l/2, 

pS(o) ^ 



/Trn 

The result of the classical walk is also in contrast with that of our model. 
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